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1999 at Langley Research Center. A copy of that briefing is attached. 

1. Unsteady Aerodynamic Modeling 

1 . 1 Fourier Functional Analysis 

The Fourier functional analysis to generate aerodynamic indicia! integrals was developed 
in Ref. 1 . The computer code to calculate the model coefficients was further improved in Ref. 2. 
However, the code still requires a user's occasional intervention to obtain reasonable results. To 
make the code more user-independent, the following improvements have been made: 

(1) The best values of the initial input variables for C L -, C D - and C m -models are now chosen 
internally in the code. 

(2) A systematic mechanism is established to allow a user to restart the calculation for a better 
solution. 

The resulting code has been applied to the F- 1 8 and F- 1 6XL data in the analysis of random 
gust response. As a result, a paper in random gust response was written and has been accepted 
for publication in the AIAA Journal (Ref. 3). 

1.2 Generalized Dynamic Aerodynamic Models 

The aerodynamic model based on indicial integrals is valid to represent the aerodynamics 
in unsteady flight conditions. However, this model may require too much computing time for real 
time simulation and does not clearly reveal the flight conditions for dynamic instability. As an 
alternative, a generalized dynamic aerodynamic coefficient model, which is frequency- and 
amplitude-dependent, is proposed. The results were published in Ref. 4. 

In addition, the developed models were also applied to the prediction of aircraft-pilot 
coupling (Ref. 5). 

1 .3 Fuzzy Logic Modeling 

The modeling methods described in Sections 1 . 1 and 1 .2 are applicable when only one 
motion variable, such as the angle of attack, appears in the model. To generalize the aerodynamic 
modeling methods to more motion variables, we developed a fuzzy logic modeling algorithm that 
could handle complex modeling problems, such as one involving the angle of attack, sideslip, 
control deflections, etc. Using this method, we developed a longitudinal unsteady aerodynamic 
model for an F-16XL configuration by combining more than 200 different data sets into one 
model. The results were published in Ref. 6. 

The same modeling technique was also applied to the rolling and yawing oscillation data 
for the F-16XL configuration. It was also possible to incorporate the rotaiy balance data into the 
present model. The rotary balance data are regarded as steady flow data and had to be represented 
as a separate aerodynamic model in the past. The results were published in Ref. 7. 
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2. Flight Dynamic Simulation 

The established unsteady Aerodynamic models for the F-16XL have been successfully 
incorporated into a conventional flight simulation code for the purpose of investigating the 
unsteady aerodynamic effects on flight dynamics. The ultimate purpose is to determine whether 
using the unsteady aerodynamics models is more realistic in design for flight safety and in pilot 
training in flight simulators. The results will be presented in August 2000 (Ref. 8). 
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Motivation and Objectives 
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Unsteady Aerodynamic Phenomena 



• Classical aerodynamic lag does not provide sufficient explanation to natural flight phenomena. 
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Classical unsteady effects - circulatory flow 



•Aerodynamic forces and moments depend strongly on the location of vortex and are therefore a function of time. 




Edge-separated vortices: 
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Overview of Unsteady Aerodynamics Modeling Techniques 
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• Quasi-Steady Aerodynamics 

• Unsteady aerodynamic effect is represented through dynamic stability parameters which are 
obtained in small-amplitude forced oscillation tests at a pre-determined frequency at a given angle of 


• Effects of "spilled vortex" and "wake vortex sheet". 

• Virtual mass effect. 

• Multiple vortices and synergistic effects. 
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• In applications, the generated indicial functions are defined as functions of a and a and interpolation 
is made to choose a particular indicial function to be used. 
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Nonlinear Indicial Function Model 
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• Nonlinear Indicial Integrals through Fourier Functional Analysis 

Refs. AIAA Journal Sept. 1992, p.2259 ° 

AIAA Paper 93-3626 

• Test data from large-amplitude forced oscillation testing in pitch with several 
oscillation frequencies are analyzed through Fourier Functional Analysis. 
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where PD’s are Pade approximants of order 2 and are defined as: 



• The effective frequency can be obtained analytically: 
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Generalized Dynamic Aerodynamic Coefficient Model 


Fuzzy Set : Allows a degree of belonging, such as 0.5 or 0.8. 
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Fuzzy Logic : 

Models the uncertainty of the definition of the event itself 
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Example : Patients suffering from hepatitis show in 60% of all cases high fever, in 45% 
of all cases a yellowish colored skin, and 30% of all cases nausea. 

High fever, yellowish colored skin and nausea are all "fuzzy". 


• The final prediction of an aerodynamic coefficient is equal to the weighted average of contributions 
of all fuzzy cells. 
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Fuzzy Logic Modeling 
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• Internal functions 

Internal functions are assumed to be linear functions of input parameters as follows 

P'=y i (x l ,x 2 ,-,x T ,-,x k ) = 



• n is the total number of cells of the model. 



• A fuzzy cell is formed by taking one membership function from each variable. The rule of the i-th 
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• The unknown coefficients of the internal functions are determined with the Newton gradient- 


• The number of membership functions (g l5 g 2 ,...g r ,...g k ) is determined by having the largest value of 
multiple correlation coefficient R 2 . 
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• The parameters associated with the internal functions are identified with Newton's method to minimi 


• “predict.f ’ - to predict the results. It requires one input file: predict.inp using the same format as for “struc.f 


3 

0 

CL 

H> 

1 

0 

3 

a* 

cd 

t 

cd 

3- 

r+ 

in 

S’ 

3 

1 

H> 

M 

? 

I 

i 

I 

f 

d 


J1 p 

o 

O 

p 

o 

o 

O 

O 

o 

o 

p 

p 

^ 00 

oo 

-O 

to 

bv 

bv 

Lh 

KM 

4^ 

4- 

to 

to 

VO 

4L 

oo 

LO 

oo 

4^ 

vO 

U> 

vO 

4^ 

vo 

4x 

h O 

o 

vo 

vO 

VO 

O 

O 

VO 

to 

O 

O 

o 

O 

C- vo 

U\ 

KM 

4^ 

4* 

U) 

u> 

to 

to 

to 

i — * 

1 — » 

S P 

p 

to 

00 

U) 

p 

4L 

p 

L/v 

' — 1 

Ov 

u> 

l "i vi 

1— » 

’ l 

b 

2 


4L 

bo 

bv 

o 

vo 

H- 4 

I—, Vi 

o 

VO 

U) 

o 

VO 

vo 

u> 

u> 

■o 

4- 

5 to 

VO 

to 

o 

o 

to 


o 

VO 

-o 

LO 

to 


& Ov 

-o 

00 

oo 

oo 

VO 

vO 

vo 

00 

00 

-o 

-0 

h-k Lh 

a oo 

CO 

o 

U1 

p 

>— 1 

to 

•— * 

vO 

4^ 

vO 

to 

-0 

4^ 

•o 

th 

to 

to 

4* 

o 

to 

to 


CD vo 

o 

00 

C>h 

4^ 

Ov 

vo 

•— * 

OJ 

to 

o 

Ov 

I to 

4^ 


00 

O 

L*J 

O 

o 

4^ 

00 


4- 


o o o o o 


X 3 

§ 

CD 

3 

I' 

Vj 

o 

g 

CD 

e- 

CD 

CL 

s* 

CD 

c/q 


O 

w 

CD 


i 


i 


H- VO -° P 

- ■- oo 2 B 

C3 00 oo w ° 


to 

00 

Vi 

00 

VO 


Vi 

4* 

VO 


o o o 


o o o 


V> Ov 

Vi Vi 

Ov to 
»— 4* 

-o o 
o o 


£ 5 t 3 £ £ 


4 * 4 * -d 


to 

4^ 


to 

4> 


to to 
4* 4* 


VO 

-■o 

to 

I— l 

-o 

o 


to 

4* 


to 

Vi 

On 

4^ 

to 


vO 

to 


VO 

VO 

H- * 

VO 

2 

4> 

4^ 

Ov 


•O 

4^ 

oo 

o 

Ov 


o o 


Vi Vi Vi 

© © © 
o o o 
o o o 

o o o 


Ul Ul Ul 


o 

o 

o 


o 

o 

o 


o 

o 

o 


o 

o 

o 

o 

b 

o 

o 


o 

o 

o 

o 

b 

o 

o 


Vi Vi Vi 

b b 
o o 
o o 


o o 

b b 
o o 
o o 


8 

O 

o 

b 

o 

o 


to 

4* 


O 

O 

O 

O 

b 

o 

o 


8 

Ov 

vO 

8 


to u> 
o to 

H- Vi 

to 
O 
o 


4^ 

O 

O 


i— • 


H-* 

Vi 

t— ‘ 

to 

pv 

1—* 

t— ‘ 


p 

b 

Vi 

Ia 

LO 

Vi 

Ov 

4^ 

o 

to 

4^ 

b 

V) 

bo 

•o 

-o 

LO 

Vi 

4^ 

o 

4^ 

4^ 

-o 

o 

O 

O 

►— * 
o 

o 

O 

VO 

o 

Vi 

o 

to 

o 

o 

o 

o 

o 


o 

o 

o 

o 


4* VO 
H- O 


Ov 4* 

oo to 
to oo 
Vi vo 


•— o 

4 ^ vO 

o o 


to o 
to vo 

4* <1 
<1 to 


Vi 4 ^ 

bo to 

m t— * 

-o 00 

H- to 
vo *“ A 
Ov 4 ^ 
to b 
-O 4* 
O Vi 


LO to 
U) 

to ov 

to to 
to <-o 

s s 

Vi b 
00 Vi 

© vo 


to 

4> 


to to 

4* 4* 


to to 

4^ 4^ 


Ul Vft V/i 


o 

o 

o 


o o 
o o 
o o 


o o o 


o 

b 

4^ 

O 


O | 

a 

© <"*■ t-H 

8 i 

n> 

J’ 

g. 

t 3 

P cti 

t- 

r-K 

O R 

< 

° o 

d* 


P o 
Ov * 0 
Ov § 

00 © 

r° 6 

£ b 

^ vi 


o o o o o 


o 

o 

o 

o 

Ov 

00 

Ov 

o 

o 

o 


o 

o 

o 


o © 
o o 
o o 

o p 
to to 
Ov 4 » 
-U Ov 

VO to 

88 


o o 
o o 
o o 

o o 

8 8 

o o 

o o 

t— l 

o 

Ov 

H- > 

o 
o 


y* 6 , 

2 © 


a* 

ft 

R - “O 

I 

1 

£1 CL 

cr 

v: 

& 

CD 

g 
3 


O 

o 


o 

o 


V o 
§ O o» 

8 g“ 


CD 

d 
O 


i' 


Ov 

2 

O 

O 


o 

b 

-o 

Ov 

■o 

o 

o 


o 

© 

o 

v» f 
oo 
o 
o 


d 

_ CL 

O 2 


CO 

9 

O 

H> 


o 
cr ~- 
?T 

co yj 

3 1 & 

a « 

I? 

5 g 

s 8 

a 

& 

CD 


<J9 

CD 

CO 

s, 


g 

a 2 


CO 

CD 


I 
< 2 
s.| 

si 5 


CO 

d 

CD 

CD 

CL 

CD 

CL 


hO 

to 


Three codes have been developed for fuzzy logic modeling: 
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Example Unsteady Aerodynamics Models 
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♦ The model structure was determined to be: 


• The predicted C L is 3.79087/8.251 18 = .4594, to be compared with .4553 in the data. 
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the following results can be obtained for the membership functions: 
A,(l) = . 1977777, A 2 (l) = .5230659, A,(2) = .4986683 
A 4 (3) = 5.57500E-02, A 5 (l) = 1.0, A 6 (2) = .5 
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• Models created with sets of data from large amplitude-oscillation tests did not predict well other sets of 
data at high frequencies, because these test conditions are not represented in these models. 













Strong effect of a on directional stability 
Correlation with training data: C n and C ( 




Directional instability becomes worse by the dynamic effect 
Predicted results: C n and C f 

k = 0.201, 5. = 0° 
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• a Effect 

• The term associated with a appears as the virtual mass effect in incompressible unsteady 
aerodynamic theory, and could be important at high frequencies. 
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• The best model structure was determined to be: 



Rotary balance data from a 10% model were incorporated into the aerodynamic model: 











Model for the Rolling Moment Coefficient - correlation coefficient - 0.948 

• C ( from roll oscillation 

• A counterclockwise hysteresis loop implies positive roll damping 




C ( from yaw oscillation 
A clockwise loop means C. is positive. 






y(deg) V(deg) v(deg) 



Model for the Yawing Moment Coefficient - Correlation coefficient — 0.80 

• Hysteresis loop of C n vs i | f curves is clockwise =* > 0. 

• C n measurements did not provide good correlation 
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• In simulation, the unsteady aerodynamic effect 
is evaluated by replacing the dynamic-derivative 
terms in the quasi-steady formulation with the 
difference between the tunnel unsteady data a nd 
the tunnel static data; 


time(sec) time(sec) 



■Dynamic derivatives from the quasi-steady aerodynamic model, with a rudder step input of -2 deg. at a trimmed 
a = 21 .4 deg. 


time(sec) time(sec) 



• Dynamic derivatives from the unsteady aerodynamic model, with a rudder step input of -2 deg. at a trimmed 
a = 21.4 deg. 
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■ Unsteady data are adjusted to correspond to an aircraft c. g. of 0.233 
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Longitudinal Stability 
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Modeling of Aircraft-Pilot Coupling 
(Ref. A1AA Paper 98-4147) 


4>*|/a>=0.38 sec 


II o 

gB* 

3 a 

U\ 

H- 


H 8 

cr b 

^ " II 

>■* 2 . h- 

ii <8 ° 

' 3 d. 

O < Ci 

tog. F 
U) P 
4^. CD 

<0 ^ ^ 

H- §» f 


o» 

I "03 


tO t/3 

oo tr 


o ^ 


Xi 

V> 

CD 

s* 

3. 

H— 1 

o 

O 

Cl 

B 

4l 

0 

V# 

o 

CL 

o 

V* 

CD 

O 

cr 

a 

V* 

O 

V* 

CD 

CL 


O 

P 

o 

hO 

VJ 

H— 1 

o 

V# 

M • 
| 

o 


o 

St 

o 

a 

S> 

V! 

& 

CD 

»t 

O 

CL 



•The following flight conditions are taken from an analysis of equilibrium surfaces and bifurcation. 
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Table 1 Effect of Time Lag on Susceptibility and Severity 
of Pilot-Induced Oscillation for the Quasi-steady Model 
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Equation of motion: 

MVa(t) + ( 1/2) p V 2 SC l {t) = ~L(t,a, .& ) 
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•Nonlinear unsteady aerodynamics produced 48.7% ~ 90.1% higher possible maximum incremental load factors 
to the von Karman random gust than the linear unsteady aerodynamics for the F- 18 HARV and the F-16XL. 



• Nonlinear stability and control analysis using simulation models with unsteady aerodynamics over the whole flight 
envelope and in abnormal flight conditions. 


